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High standards of living are linked with the availability of energy. 
Therefore, there is always a requirement to modify existing power 
generating systems or to add new systems to fulfill the increasing demands 
of energy. It is becoming a core issue in the South Asian region so that at 
least current standards of living could be maintained by controlling growing 
energy rates and shortening the availability of conventional sources. Furnace 
oil is a residue obtained in the process of distillation of crude oil in the 
petroleum industry. This furnace oil is applicable for power generation 
plants allied with the petroleum industry so that residue of oil refineries 
could be utilized. In this research, a furnace oil-based power plant is 
designed in electrical transient and analysis program (ETAP) software. 
Moreover, steady-state and transient analysis of the proposed design of the 
power plant are conducted to increase practical viability of model. Load flow 
is analyzed to depict the performance of the power plant model under load 
conditions. This model is further tested under the motor starting event to 
analyze the current drawn by the motor. This event helps to determine the 


permissible values of protective equipment installed in power plants. In the 
last, economic dispatch analysis is conducted to find the relative minima of 
generation cost of furnace oil-based power plant with respect to coal and 
hydel generation 
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1. INTRODUCTION 

With the increasing electricity demand, there is always a need to go for alteration of existing power 
generating systems and to add new generating units. These additions and modifications must be technically 
adequate and economically feasible in times of growing electricity rates and load densities. Developing 
countries like Pakistan have been facing energy shortages for the last decade. Therefore, it is the need time in 
developing countries to add new generating units so that sustainability in investment and economy could be 
obtained [1], [2]. 

Furnace oil is a residue obtained as a result of the distillation process of oil refineries. This furnace 
oil is useful for the generation of electricity. Therefore, oil refineries are highly motivated in the Middle East 
where oil reserves are very massive, to install power generating units so that residue of oil refineries could be 
utilized beneficially [3]. 

Design of power plant is linked with few analyses without which power plants could not be termed as 
functional. Among them, the load flow analysis is the key to the acquire steady state of the system [4]-[16]. It is 
not possible to conduct load flow analysis by handy calculations when the power system is big. In the 
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literature, different linear and nonlinear techniques are available to depict the performance of power systems 
under steady-state conditions. These techniques are evaluated in terms of their complexity and computational 
time. Electrical transient and analysis program (ETAP) software is equipped with the Gauss-Seidel iterative 
technique and Newton Raphson technique. In large and complex power systems, Newton Raphson is found 
more effective because of less computational time and fewer iterations for convergence. Therefore, its 
convergence rate is faster. Adaptive Newton Raphson is also available but because of zeros in off-diagonal 
entries, its convergence rate is rapid but its result loses effectivity because of lost or little information [4]-[7]. 

The efficiency and reliability of a power system can only be determined by testing its performance 
under abnormal or faulty conditions. The short circuit is a condition in which current increases rapidly with a 
sudden drop of voltage at the faulty point. The power system is required to be operated under this condition. 
Equipment involved in the power system must have the ability to sustain its operation under such 
abnormalities. Protective equipment installed must have the ability to respond to such abnormalities as 
desired. In addition to this, this analysis also helps to determine the operating limits of the system, sizing of 
equipment, and rating of power system components [11], [17]-[19]. 

The phenomenon of motor starting is also an important event in the power system because of 
starting characteristics of the motor. 85 percent of motors installed in power systems are induction motors. It 
draws six to eight times larger current as compared to its rated current at the time of starting. An increase in 
current makes this event very similar to the event of a fault. Therefore, it needs to be discriminated among 
these events [20]-[23]. 

There are events in the power system which increase current for a very short interval of time. These 
events are mostly termed transients. These events decay down on their own or it needs just a restart of the 
system at a low scale. Therefore, steady-state events must be differentiated from transient events as steady- 
state events are required to be rectified [24]-[39]. 

Cost of generation is an important pillar of any power generating unit. New installations must be 
reasonably economical so that consumers and supply company's relationship could be derived in an optimum 
way. Economic dispatch is a way to determine the optimal generating cost from a particular plant while taken 
into consideration other power generating units. This analysis helps to determine the operating conditions of 
the power plant under minimum fuel cost [40], [41]. In this research, interior point optimization is applied 
because of its rapid convergence in large power systems and is an effective way of handling large power 
systems in ETAP software [40]-[50]. 

There are different types of softwares available for the simulation of models of power systems and 
power stations. In all of these softwares, ETAP software is found most promising because of its realistic 
approach and user's control of design and analysis parameters. Therefore, ETAP is associated with many 
industries for simulation and development [28], [51]-[59]. 

In this research, a furnace oil-based power plant consisting of more than 160 bus bars of high 
voltage, medium voltage, and low voltage is designed. Moreover, post-design analyses discussed in the 
coming section are performed to depict this power plant design performance. This research paper consists of 
the following sections, Section I covers the introduction and literature review. Section II consists of the 
methodology involved for the design and analysis of furnace oil-based power plants comprehensively. 
Section III contains the test system designed for simulation under different scenarios. Section IV consists of a 
discussion of the results of the simulation. The conclusion is presented in Section V. 


2. METHODOLOGY 
2.1. Data collection 

Relative data for the design of the power plant is gathered. A single line diagram should be provided 
by the vendor. If collected data is sufficient for the power system studies, then electrical design and analysis 
program (ETAP) software is used to design a model. If some ratings are not correct or it is not available, the 
site is visited to gather information. After this, if still data is not sufficient for the model which is collected 
during the site survey, then assumptions are made based on the different international standards such as IEEE 
standards, ANSI standards, IEC standards. Different analyses are conducted using ETAP to generate a 
recommendation report for the power plant. Load flow analysis (LFA) is the first and the foremost technique 
of system studies. The results of this study further help in the implementation of techniques. 


2.2. Load flow analysis 

Load flow analysis is the study of power, voltages, and angles at different nodes of bus and specify 
that whether the system is operating in normal voltages condition or emergency operating conditions, Load 
flow analysis is performed to design a new power system and modification of the existing system. The 
voltage (V), voltage angle (ô), real power (P), and reactive power (Q) can be computed using load flow. 
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Three types of buses are used in load flow which are swing, generator, and load buses. In swing bus, P and Q 
are found based on and V and ò. In the generator bus, for the input of P and V, the required output of Q and ô 
is calculated, and in the load bus, P and Q are known and V and 6 are unknown values. All the unknown 
parameters (V, ò, P, and Q) are accessible for swing bus, generator bus, and load bus with load flow analysis. 
The other electrical parameters such as current, power factor, apparent power for the whole electrical system 
can be determined based on these four parameters (V, 6, P, and Q) and impedance between the buses (Zij). 
For the calculation of load flow analysis, there are some methods i.e. Gauss-Seidel method, Newton Raphson 
method, Adaptive Newton Raphson method. In this research, the Newton Raphson method is implemented 
because its results are accurate and it requires fewer iterations. Basic model of power system is presented in 
Figure 1 in which current Ii is entering from the bus (Vi) and is flowing towards different buses (V1, V2, ..., 
Vn) via admittances (yil, yi2, ..., yin). 


Figure 1. Basic power system 


2.2.1. Newton raphson 

Newton's method, named after Isaac Newton and Joseph Raphson, is a method for finding 
successively better approximations to the roots of a real-valued function. The approach to Newton-Raphson 
load flow is similar to that of solving a system of nonlinear equations using the Newton-Raphson 
method [4]-[7]. It is an iterative method that approximates a set of non-linear simultaneous equations to a set 
of linear simultaneous equations using Taylor's series expansion. Referring to Figure 1, power flow 
equations are formulated in polar form for the n bus system in terms of bus admittance matrix Y as: 


l= diet YV; a) 
Where ith and jth buses are donated by i and j. Expressing it in polar form as: 

I; = Xj-1 [Yy lI | < Oy + ô; (2) 
The current is expressed in terms of the active and the reactive power at the bus Z as: 

jez (3) 
Substituting for Ji from (2) in (3) as: 

P, — JQ; = [Vil € —8 dina IVil[Vy] < 6i; + 6; (4) 


Separating the real and imaginary parts as: 
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Qi = — YL, VViVjYij|sin (Bij 8j — di) ©) 


After expanding (5) and (6) in Taylor's series with initial condition and neglecting higher-order terms as: 


OP. OP. OP. OP. 

d ger ee 2 — ? lr Ad, AP; 

86; ðn Val EIN, : 

OP OPn ôP OP, : : 

35, | an Ov “lal |] A0, | |AF, (7) 
92 80; Qz ^ 90» ||A|V;| AQ; 

862 ðn Va} 8|V4] : : 

9Qn Qn In 8Qn AIV, | AQ 

08; ðn ôV "7 Val n n 


The Jacobian matrix gives the literalized relationship between small changes in Ad; and A|V;| voltage angle 
and magnitude with the small changes in real and reactive power AP; and AQ; as: 


hol= b nlar ® 


The diagonal and the off-diagonal elements of J1 are: 


OP; 


a; = XjaVivjYül cos(6ij + 6j — ài) (9) 
J 

22: = Y? |VifYij| sin (ij + 8j — ôi) (10) 
J 


Similarly, the diagonal and off-diagonal elements of J2, J3, and J4 can be evaluated. The terms ap," and 
0Q;* are the difference between the scheduled and calculated values, known as the power residuals given as: 


APF = P” — PE (11) 


AQF = QF” — Qk (12) 


Using the values of the power residuals and the Jacobian matrices, Ad; and A|V;| are calculated from (7) to 
complete the particular iteration and the new values are calculated as shown below which are used as the 
initial values for the next iteration. 


6/41 = eë — ASF (13) 
[Vi] = [VF] — [AVE | (14) 


2.3. Short circuit analysis 

Short-circuit studies are performed to determine the peak magnitude of the prospective currents 
flowing in the power system at different time intervals after a fault occurs. The magnitude of the currents 
flowing through the power system after a fault varies with time until they reach a steady-state condition. This 
behavior occurs due to system characteristics and dynamics. During this time, the protective system is called 
on to detect, interrupt, and isolate this faulty equipment. The duty enforced on this equipment is dependent 
upon the magnitude of the current, which is dependent on the time from fault inception. This is done for 
different types of faults (three-phase, phase-to-phase, double-phase-to-ground, and phase-to-ground) at 
dissimilar locations throughout the system. The information is used to select the size of the equipment like 
fuses, breakers, and switchgear ratings in addition to setting protective relays. 

Motor starting analysis, Motor starting analysis is used to calculate motor acceleration time and the 
voltage impact on motor starting in electrical power systems. Advanced plotting and time-varying graphical 
displays enable engineers to quickly evaluate results and to make decisions. In an induction motor, when 
supply is given to the stator windings, the rotating magnetic field flux, and the produced flux in the rotor 
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windings due to the back emf, cause the motor torque to increase resulting in high rotor current. During the 
time between the application of electric supply to the motor and the actual acceleration of the motor to its full 
speed, a large amount of current is drawn by the stator from the supply. This starting current is about 5 to 6 
times more than the full load current. This time duration can be for a few milliseconds or longer. This causes 
the electrical equipment to damage because of the increased voltage drop in electrical systems due to the flow 
of larger currents across the cable. For this reason, a definite method of starting the motor is needed. 


2.4. Transient analysis 

The capacity of a power system, containing at least two synchronous machines, to keep on working 
after a change happened on the system is a check of its stability. The stability issue takes two structures: 
steady-state and transient. Steady-state security might be characterized as the capacity of a power system to 
keep up synchronism between machines inside the system following generally moderate load changes. 
Transient stability is the virtue of the system to stay in synchronism under transient conditions, i.e., abnormal 
and switching operation, etc. In a modern power system, stability may include the power organization system 
and at least one in-plant generator or synchronous machine. Following are the possibilities i.e., for example, 
load dismissal, unexpected loss of a generator or utility tie, starting of motors or fault (and their time 
interval), directly affect system stability and efficiency. Load-shedding schemes and basic flaw clearing 
times can be resolved to choose the best possible settings for protection relays. These kinds of studies are 
likely the absolute most complex ones done on a power system. The beginning stage for any transient 
stability simulation is the pre-disturbance control power flow arrangement. The system is constantly thought 
to be in a consistent state when the aggravation happens. Utilizing the pre-disturbance power flow solution, 
all the underlying estimations of the state factors are determined after these are acquired, the disturbance is 
simulated. The simulation of the disturbance causes a conflict between the mechanical power contribution to 
the generators and the electrical power yield of the generators. The equilibrium is disturbed. The state factors 
calculated by the different differential conditions combined with mathematical conditions change their 
values. A transient stability study has been done by introducing a three-phase and L-G fault and its clearance 
on different buses. The system is analyzed under the transient produced by these abnormal conditions. The 
system generator's reaction at beginning of this issue and clearance of this disturbance has been recorded. 


2.5. Economic dispatch 

Due to the increasing price of fuel oil, wages, and salaries, the electrical energy cost also increases, 
so one should keep in mind the cost factor while operating or building a new power system, to meet the 
growing load demand. In the operation of the power system, the contribution from each load and each unit 
within a plant must be such that the cost of electrical energy produced is minimum. Interior point 
optimization is the technique through economic dispatch is conducted in ETAP. 

The reduction of generation cost is called economic dispatch as shown in Figure 2. In the 
generation, cost includes the operational cost, maintenance, labor, and supply cost. Mostly, the generation 
cost is present as quadratic of input fuel and power generation cost in economic dispatch and other 
constraints are voltage limits, real and reactive power limits, and power balance condition [40]-[45]. To 
understand the economic dispatch, a problem is presented in the mathematical description as: 


Min XP (e; + BiPoi + Yi Pao) 
Subject to: 
XP) = Pp +P, (15) 
PRY < Poi < PRA 
Where a, f, and y are the parameters that are dependent on the system and they affect real generation power 
cost. The total number of generators is represented byNG . P; is the real power generation at bus i and Pp and 


P, are total power demand and power losses in the power transmission network respectively. The quadratic 
interior-point program converts economic dispatch problem in the following format [46]-[50]: 


Min (AT E - BTx + : xT Qx) 
Subject to 


pmin < Ax < pmax (16) 
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F(x)=0 
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Figure 2. Proposed flow chart of inter-point optimization for economic dispatch 


3. TEST SYSTEM 

The test system is composed of three main sections which are described below. It contains nine 
furnace oil-based generators and one combined cycle-based generator. Further detail is followed by each 
section. 


3.1. High voltage side 
The high voltage side is composed of two main buses which are directly fed from two utility buses 
as given in Figure 3. These buses have a 132 KV nominal voltage level. 


3.2. Medium voltage side 

Voltages are stepped down from 132 to 15 KV for connecting it to generators as shown in Figure 4. 
The medium voltage side is composed of three main generator buses. Gen bus-1and Gen bus-2 are fed from 
generators 1,2,3,4 and generators 5,6,7,8 respectively. Each generator has a rated capacity of 17500 KW. 
Every generator bus is linked with each other using a bus coupler to maintain the supply in case of any fault. 
Gen bus-3 is connected directly to Generators 9 and 10. Gen 9 and 10 have rated capacity of 17500 and 
10750 kW respectively. Now again, voltage levels are stepped down to 0.415 KV for supplying to the LV side. 
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3.3. Low voltage side 


Low voltage is modeled by using a composite network for avoiding the complexity in the case of a 
large-scale power network. As explained before, low voltage side buses have a nominal voltage level of 
0.415 kV. This section is consisting of three feeders as shown in Figure 5, Figure 6 and Figure 7. 

Each feeder is further composed of sub-buses. These buses are attached to the auxiliary loads 
required for operating the whole power unit. Auxiliary loads have a long list of loads. So, the loads are 


lumped on LV buses depending upon their rated kW. All the loads which are less than 30 kW are lumped 
while modeling on ETAP. 
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Figure 4. A sectional view of medium voltage side 
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Figure 6. A sectional view of low voltage side feeder 2 
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Figure 7. A sectional view of low voltage feeder 3 


4. SIMULATION RESULTS 
4.1. Load flow analysis 


A load flow calculation determines the state of the power system for a given load and the 
distribution of generation. It represents a steady-state condition as if that condition has been held fixed for 
some time. In reality, power flows and bus voltages fluctuate constantly by small amounts because loads 
change constantly as lights, motors, cooling equipment, and other loads are turned on and off. However, 
these small fluctuations can be ignored in calculating the steady-state effects on system equipment. Load 
flow study is performed on the basis and recommendations of the following international standards: IEEE 
STD 399-1997. (IEEE recommended practice for industrial and commercial power system analysis) and 
IEEE STD 141-1993 (IEEE recommended practice for electric power distribution for industrial plants). After 
performing load flow, voltages at each bus are checked which should be under safe level as shown in Table 1. 


Table 1. Percentage loading on each bus 


Bus Nominal Voltage (kV)Voltage (%)Nominal Current (A)Loading (A)% Loading 
MAIN-1 132 100 2000 402.3 20.1 
MAIN-2 132 100 2000 342.8 17.1 
G-BUSI 15 105.86 3600 2532 70.3 
G-BUS2 15 105.85 3600 2531 70.3 
G-BUS3 15 102.96 3600 1085 30.1 
ENG-WISE-CONNSUMERI 0.4 104.79 5000 1567 31.3 
ENG-WISE-CONNSUMER2 0.4 103.23 5000 1429 28.6 
ENG-WISE-CONNSUMER3 0.4 97.79 5000 1578 31.6 


Figure 8 explains the graph between the nominal and actual voltage on different buses. As it is seen 
that loading at mainland main 2 buses are within the limits of 95-105% so it depicts that all the buses are 
sized properly. Figure 9 graph between the maximum withstands capacity of buses for short circuit fault and 
actual fault current flowing through the buses. All the buses have a fault current that is lesser than the 
maximum withstand capacity of the buses. For example; G bus! and G bus 2 have fault current of 2.531 kA 
and 1.0835 kA respectively and their maximum withstands capacity are 3.6 kA respectively. 

Figure 10 is the graph between the rated power and actual active power of each generator at a 
constant power factor. As it is seen from the graph that Gen 7 and steam turbine generator have rated active 
power of 17.5 MW and 10.75 MW respectively where their actual active power is 16 MW and 9 MW 
respectively. One thing that is noted is that the required MW is attainable from the generators at a power 


factor of 0.85. Table 2 covers the details of rated active power (KW) and rated reactive power (KVAR) of 
each generator. 
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Figure 8. The graph between the nominal and actual voltage on different buses 


OLoading current (KA) 


ElNominal current (KA) 


sist [ 
s Ml 


eri L —— 
s lll! 


ist Co 
s Il! T! 


ssor [— 
se mumm 


west O 
se mmm 


ust O 
se mumu 


szreo O 
c o 


zoro C 
z mmm 


“4 
Vb 
+ 


^4 
Y). 
4 


Figure 9. Graph between nominal and loading current of different buses 
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Table 2. Rated KW and KVAR of each generator 


Generators Nominal Voltage Rating (kW) KW KVAR % PF %Generation 
(kV) 
Gl 15 17500 16000 10845 85 100 
G2 15 17500 16000 10845 85 100 
G3 15 17500 16000 10845 85 100 
G4 15 17500 16000 10845 85 100 
G5 15 17500 16000 10845 85 100 
G6 15 17500 16000 10845 85 100 
G7 15 17500 16000 10845 85 100 
G8 15 17500 16000 10845 85 100 
G9 15 17500 16000 10845 85 100 
Steam Turbine 15 10750 9000 6663 85 75.8 
Generator 


4.2. Short circuit analysis 

The graph depicts the relationship between short circuit current, three-phase fault current, the line to 
ground, line to line current with respect to the maximum withstand capacity of each bus. As seen from 
Figure 11 that at Gbus-3, three-phase device duty, a line to ground, and line to line faults have the values of 
35.1, 30.74, and 30.58 KA respectively. The maximum withstand capacity of G bus-3 is 40KA as shown in 
Figure 11. So, from the observations, it is concluded that the different fault values exist at different buses. 
Based on these values, the protection scheme for the tested system is designed. 
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Figure 11. Graph of short circuit current, three-phase fault current, the line to ground, line to line current 


4.3. Motor starting analysis 

The test system is analyzed by starting the biggest motor of the system and its effects on all the 
buses are observed and are presented in Table 3. This motor is connected with the black start unit bus. After 
performing motor starting analysis, all the profiles are checked with respect to the time which is within the 
safe limits. Voltage dip is seen from the graph on the black start unit bus. The graph shown in Figure 12 
explains the variation in voltage at different buses with respect to time. 

It is seen from the graph that at t=1sec, the motor starts and there is seen a voltage dip at Main 1 and 
Main 2 which is 131.986 KV when the nominal voltages on both buses are 132 KV as shown in Figure 12. 
There is a voltage dip 0.2 KV at the main bus. It is also seen from the graph that after few seconds, voltage 
again retains its normal value which is 132 KV. It is concluded that voltage drop is within the marginal limits 
which are followed internationally for the motor starting purpose. It is noted that motor starting analysis is 
done with the respect to 50 kW motor starting at black start unit bus. Figure 13 shows the ETAP generated 
graph of bus voltage with respect to time. At the black start unit bus, there is a voltage dip at 1 second when 
the motor starts and after some time voltage again retains to its level. All the other buses have no dip because 
the corresponding motor is connected to the black start unit only. 
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Figure 12. Graph of variation in voltage at different buses with respect to time 


Table 3. Voltage profiles at various intervals analyzed during motor starting analysis 
Nominal Voltage (kV) Voltage profile (%) at various time intervals (Seconds) 


Bus @ t = 0+ (Pre-start)@ t = 1(At start-up) t23(Post-start) 

MAIN 1 132 100 99.99 100 

MAIN 2 132 100 99.99 100 
Gbusl 15 105.86 105.81 105.86 
Gbus2 15 105.86 105.81 105.86 
Gbus3 15 102.96 102.96 102.96 
Eng-Wise-Connsumerl 0.4 104.79 104.75 104.79 
Eng-Wise-Connsumer2 0.4 103.47 103.41 103.47 

Eng-Wise-Connsumer3 0.4 97.79 97.71 97.79 
Black start unit 0.4 103.47 101.28 103.23 
Bus Voltage 
— BLACK START UNIT — ENGINE WISE CONSUMERI— ENGINE WISE CONSUMER 3— ENGINEWISE_CONSUMERS2- Gbuzl — Gbus2 — Gouss 
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Figure 13. ETAP generated a graph of bus voltage with respect to time 
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4.4. Transient stability analysis 


Transient analysis is performed under two cases. kW, kVAR, Terminal current, Speed, and other 
parameters of generators (presented in Figure 14) are analyzed. 
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Figure 14. Graph of generated MW of different generators with respect to time 


a. Case Sc-GB 

In this case, a three-phase fault is introduced on generator bus 2 and is cleared after some interval. 
Figure 14 shows the graph of generated MW of different generators with respect to time. It is seen that at 
Gen1 and Gen2. At t=0 and t=3, the active power is 161.32 MW which remains constant when a transient is 
induced on the generator at bus 2. Figure 15 shows the Graph of reactive power with respect to time. It is 
easily cleared from the graph that at reactive power of generators 7 and 9 have transient at t=3sec. After some 
time, the transient is again normalized to the level before the transient. It is seen that the transient time is in 
the constraint which is internationally followed. Figure 16 is the graph of terminal current with the variation 
in time. It is seen from the graph at t=3sec generator 9 have a large terminal current with respect to the other 
generator such as generator 2,3and 4. Gen 9 terminal current is 3000A and generators 2, 3, and 4 have 
terminal currents of 500A respectively. When a three-phase fault is introduced on G-Bus2, then the system 
will operate in a transient state as seen in Figure 17. Voltage dip is seen at the G-buses. This abnormal 
condition is settled in some period of time and the operation of the generator also becomes normal. 

In transient conditions, all other parameters are disturbed but there is a significant effect of transient 
conditions on generators MW. As seen in Figure 18, MW of Gen-9 is raised to a very high value of 50 MW 
as compared to the other generators. It is noticed that the effect of transient is more prominent on Gen -9 with 
respect to the other generators but the value will remain in considerable range. A protection system will 
correspond to these types of transient very efficiently. This attribute will increase the reliability of the system. 
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Figure 15. This ETAP-generated graph is between MVAR and time 
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Generator Terminal Current 


Figure 16. This is the plot of the generator terminal current with respect to time 


Bus Voltage per Hz 


Generator Electrical Power 


Figure 18. This figure shows the generator’s MW with respect to time 


b. Case Sc-GB2 

In this case, a line to ground fault is introduced at generator bus2 and removed after some interval. 
After some observation and working on the tested system, a comparison is established with respect to the 
time of transient case-2. In transient case-2, a line to ground fault is induced on G-bus 2. As shown in Figure 
19, at t=3.5 seconds, Gen 9 has an active power rating of 2.034 MW which is a huge variation. This behavior 
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illustrates the effect of transients very clearly. Figure 20 shows the transient stability in terms of rated active 
power (MW). Stability is achieved in 30 seconds. In 20 seconds, transient stability is achieved in term of 
rated reactive power (MVAR) as shown in Figure 21. Rate of change of active power and reactive power 
ultimately affect the current and voltage. Therefore, under transient stability analysis, current and voltage 
study are also done and is presented in Figure 22 and Figure 23. 
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Figure 20. This is the plot of MW with respect to time in transient case-2 
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Figure 21. This is the plot of MVAR and time in transient case 2 
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Figure 22. This is the graph of generator terminal current and time in transient case 2 
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Figure 23. This is the graph between bus voltages and time in transient case 2 


4.5. Economic dispatch 

In the end economic dispatch of the tested power system is done using the optimal load flow module 
of ETAP as shown in Figure 24. Optimal load flow ensures the usage of maximum power with the variable 
amount of constraint which could be varied according to the different priorities and scenarios. Economic 
dispatch study is performed in a larger network where a furnace-based power plant was generating electricity 
along with a coal-fired and a hydropower plant. It appears from the study that it is economically optimum to 
utilize a furnace oil-based power plant at 56.7% of its total capacity. This also emphasizes the need to utilize 
furnace oil power plants even though cheaper fuel power plants are present in the system. 
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Figure 24. The graph between the cost of different fuels wrt to the MW at simple and optimum load flow 
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5. CONCLUSION 

Growing demands of energy can only be fulfilled with the addition of new power generating units 
and modification of existing power generating systems. In addition to this, there is a need to look up new 
energy sources to meet the energy requirement in an optimum way. In this comprehensive research, a furnace 
oil-based power plant is developed and analyzed under different pre-design and post-design analyses so that 
recommendations could be prepared under international standards. Therefore, load flow analysis, short circuit 
analysis, motor starting analysis, and transient analysis are conducted to depict the practical feasibility of 
equipment involved in the design of furnace oil-based power stations. Because of conducting simulations 
under IEEE standards, a very practical approach is followed. In addition to this, ETAP also adds fuel to this 
research because of its practical viability. Economic dispatch is also performed to find out the optimum 
rating capacity of a power plant with minimum fuel cost in an interconnected system. The compatibility of 
this research with 
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